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ABSTRACT 



We report the observation of a very high energy 7-ray source, whose position 
is coincident with HESS J1841-055. This source has been observed for 4.5 years 
by the ARGO-YBJ experiment from November 2007 to July 2012. Its emission is 
detected with a statistical significance of 5.3 standard deviations. Parameterizing 
the source shape with a two-dimensional Gaussian function we estimate an ex- 
tension a = (0.40t[J;!) , consistent with the HESS measurement. The observed 
energy spectrum is dN/dE = (9.0±1.6) x 10" 13 (£ /5TeVy 2 - 32±0 - 23 photons cm" 2 
s _1 TeV -1 , in the energy range 0.9-50 TeV. The integral 7-ray flux above 1 TeV 
is 1.3±0.4 Crab units, which is 3.2±1.0 times the flux derived by HESS. The dif- 
ferences in the flux determination between HESS and ARGO-YBJ, and possible 
counterparts at other wavelengths are discussed. 

Subject headings: gamma rays: general 



1. Introduction 



Very High Energy (VHE) 7-ray astronomy opened a new window to explore the extreme 
non-thermal phenomena in the universe. VHE 7-rays are tracers of non-thermal particle ac- 
celeration and are used to probe the conditions and the underlying astrophysical processes 
inside their sources. In the past decade, great progresses have been made in the field of VHE 
7-ray astronomy. More than one hundred VHE 7-ray emitters have been detected, belong- 
ing to several categories: Active Galactic Nuclei (AGNs), Pulsar Wind Nebulae (PWNs), 
SuperNova Remnants (SNRs), X-ray Binaries (XBs), and starburst galaxies. However, there 
is a fraction of VHE sources still unidentified because do not appear to have obvious coun- 
terparts at other wavelengths. This kind of sources may constitute a new class of objects 
with different emission properties. 

HESS J1841— 055 is an enigmatic unidentifie d VHE 7-ray source dis covered by the HESS 
collaboration during the Galactic plane survey ( lAharonian et al.l 120081 ). Its image shows a 
high extension, the measured axes for an elongated two-dimensional Gaussian shape being 
0.41° ± 0.04° (major) and 0.25° ± 0.02° (minor). HESS J1841-055, therefore, is one of the 
most extended sources in the VHE 7-ray band. The spectrum is best fitted by a simple 
power law with photon index a = —2.41 ± 0.08 in the energy range from 0.54 TeV to 80 
TeV. The integral flux is 9.1 x 10~ 12 photo ns cm -2 s _1 at energies above 1 TeV, about 40.3% 
of the Crab unit (lAharonian et al.ll2006ai ). 
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To date, no obvious counterpart has been found at other wavelengths. The wide VHE 
7-ray morphology sugg ests that HESS J1841— 055 may be the blend of multiple sources. 



Aharonian et al.l (120081 ) found four candidates which could be responsible for at least part of 



the entire VHE 7-ray emission: the two pulsars PSR J1841-0524 and PSR J1838-0549, the 
diffuse source G26.6— 0.1, which is a candidate SNR basing on its ASCA spectru m, and finally 



the hi gh-mass XB AX J1841. 0—0536. Basing on a striking spatial correlation. ISguera et al. 



(120091 ) propose that the Supergiant Fast X-ray Transient (SFXT) AX J1841.0-0536 could 
be responsible for at least a fraction of the VHE 7-ray emission from HESS J1841— 055, 
thus being the prototype of a new class of Galactic tra nsient MeV/TeV em i tters. Using 
7-rays with energies > 100 GeV detected by ,Fer mi-L AT. iNeronov fc Semikoz I (120101 ) found 
an event cluster adjacent to HESS J1841— 055, and a more extended event cluster at the 
opposite side. This may be an evidence that HESS J1841— 055 is co mposed of at least two 
different components. On the other hand, iNeronov &: Semikoz I (120121 ) suggest an association 
only with PSR J 1841— 0524, which is situated in the center of the extended source. 

The ARGO-YBJ experiment is an air shower array with large field of view (FOV) 
which continuously monitor the northern sky. The emission from the Crab Nebula has been 
detected with a statistical significance of 17 standard deviations (s.d.) at energies around 
1 TeV. With such a sensitivity, other 4 known VHE 7-ray sources ha ve been detected wit h 
significance greater than 5 a: Mrk 421 (IBartoli et alJboilah. Mrk 50 1 (jBartoli et al.ll2012al ). 
and the two extended sources MGRO J2031+41 jBartoli et al.ll2012bl ) and MGRO J1908+06 
(IBartoli et al.ll2012a ). It should be pointed out that the fluxes of the two extended sources 
measured by the EAS arrays Milagro and ARGO-YBJ are much higher than that determined 
by the Cherenkov arrays, showing that there are some systematic di fferences between the two 
observation techniques for extended sources JBartoli et alJl2012bl fl I Abdo et al.ll2012l ). Since 
also HESS J1841— 055 is an extended source, its study would benefit from an observation 
using EAS arrays. HESS J1841-055 is observed by ARGO-YBJ, at the edge of its FOV, 4.8 
hours per day with zenith angle less than 50°, culminating at 35.7°. This work presents the 
observation results for HESS J1841— 055 with the ARGO-YBJ experiment. 



2. The ARGO-YBJ experiment 



The ARGO-YBJ experiment is a full coverage extensive air shower array resulting from a 
collaboration between Chinese and Italian institutions and is designed for VHE 7-ray astron- 
omy and cosmic ray observations. The detector is operating at the Yangbajing International 
Cosmic Ray Observatory (Ti bet, P.R. China) , at an altitude of 4300 m a.s.l.. The detec- 
tor, extensively described in (jAielli et al. 2006 . 2009c ). consists of a single layer of Resistive 
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Plate Chambers (RPCs, 2.8 m xl.25 m), equipped with 10 logical pixels (called "pads", 
55.6 cm x 61.8 cm) used for triggering and timing purposes. 130 clusters (each composed 
by 12 RPCs) are installed to form the central carpet of 74 m x 78 m with an active area of 
~93%, surrounded by 23 additional clusters ("guard ring"). The total area of the array is 
110 m x 100 m. The arrival time of the particles is measure d by time-to-digital converters 
(TDCs) with a resolution of about 1.8 ns (jAielli et al.ll2009d ). To calibrate the 18,360 TDC 
channels, a software method has been developed using cosmic ray showers ( He et al. 12007 ). 



The c alibration precision is 0.4 ns and the procedure is applied every month (jAielli et al. 
2009ah . 



The central 130 clusters started taking data in July 2006, while the complete ARGO- 
YBJ detector including the "guard ring" collected data since November 2007. The RPC 
carpet is connected to two ind ependent data acq uisition systems, corresponding to the shower 
and scaler operation modes ( Aielli et al.l 120081 ). In the current work, only data from the 
shower mode are used. In shower mode, the ARGO-YBJ detector is operated by requiring 
at least 20 fired pads (N pa d) within 420 ns on the entire carpet detector. The trigger rate is 
3.5 kHz with a dead time of 4% and the average duty-cycle is higher than 86%. 

The high granularity of the apparatus allows a complete and detailed space-time three- 
dimensional reconstruction of the shower profile and therefore of the incident direction of the 
primary particle. Through the analysis of the position, size and shape of the reconstructed 
Moon and Sun shadows in the cosmic ray flux, the angular resolution, p ointing accuracy and 
stability of the A RGO-YBJ detector array have been thoroughly tested (jBartoli et al.ll2011bl : 
Aielli et al.ll201ll ). The Point Spread Function (PSF) is quantified using a parameter ip m as 
the opening angle containing 71.5% of the events. For cos mic ray-induced air s howers ipjo 
is 2.8° for N pad ~ 20, while becomes 0.47° for N pad > 1000 JBartoli et al.l boilal lbh . in good 
agreement with Monte Carlo predictions. The simulations show that the angular resolution 
for 7-induced showers is 30—40% smaller. The effective area of the detector for 7-induced 
showers depends on the 7-ray energy and incident zenith angle, e.g., it is about 10 m 2 at 
100 GeV and > 10,000 m 2 above 1 TeV for a zenith angle of 20° JAielli et al.ll2009bh . 



3. Data analysis 



The data set used in this analysis refers to the period from November 2007 to July 2012. 
The total effective observation time is 1492.6 days. For the analysis presented in this paper, 
only events with a zenith angle less than 50° are used, and the data set is divided into six 
groups according to N pa d. To achieve a good angular resolution, the event selections used in 
(jBartoli et al.ll2011al ) are applied here. The total number of events after filtering used in this 
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work is 2.42 x 10 11 . The opening angles ip-ro for events with N pa d > 60 and N pa d > 100 are 
1.36° and 0.98°, respectively. For the data set in each group, the whole sky map in celestial 
coordinates (right ascension and declination) is divided into a grid of 0.1° x 0.1° bins and filled 
with detec ted events according to their reconstructed arrival direction. The "direct integral 
method" ( jFleysher et al.l 12004 ) is adopted to estimate the cosmic-ray background and to 
ext ract the excess of 7- induced showers from each bin. The correction procedure described 
in ( iBartoli et al.ll2011al ) has been applied to remove the effect of cosmic ray anisotropy on a 
scale of 11° x 11°. A Gaussian smoothing method is used to take into account the PSF of 
the ARGO-YBJ detector. That is, the events in a circular area centered on the bin with an 
angular radius of 1. 3^70, are s ummed after weighting with the Gaussian-shaped PSF. The 
Li-Ma method (ILi &: Malll983l ) is used to estimate the significance of the excess in each bin. 



With this procedure, the northern sky has been surveyed (jCao &: Chen! l201ll ). The 
significance of the excess observed from the direction of the Crab Nebula is 17 s.d., indicating 
that the cumulative 5 s.d. sensitivity of ARGO-YBJ has reached 0.3 Crab unit for point 
sources. The sensitivity is dependent on the declina tion of the source, being degraded by 
a factor 3.5 at the declination of HESS J1841— 055 (jCao fc Chen! 120111 ). For an extended 
source with a symmetrical two-dimensional Gaussian shape with a = 0.40°, the sensitivity 
is degraded by 15%. Therefore, a simple estimation indicates that the flux from HESS 
J1841— 055 should be about 1.2 Crab units in order to be detected by ARGO-YBJ with 
5 s.d.. The required flux slightly varies if the spectrum is different from that of the Crab 
Nebula. 



4. Results 

The significance map around HESS J1841— 055, as observed by ARGO-YBJ using events 
with Npad > 60, is s hown in Figure 1. For comparison, the twelve sources in the second 



Fermi-LAT catalog (INolan et al.l |2012| ) around HESS J1841— 055 are also marked in the 
figure. Weak excesses are observed along the Galactic plane, indicating a diffuse 7-ray 
emission. A n analysis of the diffuse 7-ray emission using ARGO-YBJ data can be found in 



(IMall201lh . The highest significance is 5.3 s.d. at 0=18^39™ and <5=-6°3' (J2000), which is 
displaced 0.7° from the center of HESS J1841— 055. To estimate the statistical error of the 
position, the data are sampled 20,000 times and the statistical errors in both directions are 
about 0.45°. However, most of the excesses overlap the extended region of HESS J1841— 055 
and its gravity center (a=18 h 40 m ±12 m and 5=-5°52'±13'), obtained using all the pixels with 
significance greater than 3 s.d. within 3° x3° around HESS J1841— 055, is 0.4° off the center of 
HESS J1841— 055. These displacements may be caused by different concurring effects beside 
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fluctuation. (1) Complex morphology. According to the HESS result, HESS J1841— 055 
possibly has two or three peaks and the positions of the two largest ones are both 0.44° 
off the center. (2) The systematic pointing error of ARGO-YBJ is 0.2°, slightly increasing 
at the boundary of the ARGO-YBJ FOV. (3) The contribution of the nearby VHE source 
HESS J1837— 069, partially containing its emission. Therefore, the signal position observed 
by ARGO-YBJ largely overlaps HESS J1841-055. 

The intrinsic extension of HESS J1841— 055 is determined by fitting the distribution of 
9 2 for the events exceeding the background as shown in Figure 2, where 9 is the angular 
distance of each event to the position of HESS J1841— 055. To achieve a good angular 
resolution, only events with N pa d > 100 are used in this fit. In order to fit the data, a set of 
7-rays is generated taking into account the Spectral Energy Distribution (SED), the intrinsic 
source extension, and the detector PSF. The extension is estimated by minimizing the \ 2 
between data and generated events, from 0° to 1° with steps of 0.1°. Assuming a spectral 
index —2.3, the intrinsic extension is determined to be a ex t = (O^O^q 3 , 2 ) . It is found that 
the dependence on the SED is negligible within the uncertainties. This result is consistent 
with the estimation by the HESS colla boration, i.e., 0.41° ± .04° and 0.25° ± 0.02° along 



the major and minor axes, respectively (lAharonian et al.ll2008l ). 



Assuming an intrinsic extension a ex t = 0.40°, we estimate the spectrum of HESS 
J1841— 055 using th e ARGO-YBJ data with the conventional fitting method described in 



( IBartoli et al.ll2011al ). In this procedure, the expectation function is generated by sampling 
events in the energy range from 10 GeV to 100 TeV and taking into account the detailed 
ARGO-YBJ detector response, assuming a power law with its spectral index as a parameter. 
We define five N pad intervals: 60-59,100-199, 200-499, 500-999, and >1000. The best fit 
to the SED and the corresponding 1 a error region are shown in Figure 3. The differential 
flux (TeV -1 cm~ 2 s _1 ) in the energy range from 0.9 TeV to 50 TeV is 

^ = (9.0 ± 1.6) x 10- 13 (£/5 TeV)- 2 - 32±0 - 23 . (1) 
dE 

The median energies of the five N pa d intervals are 2.3, 3.5, 7.1, 14 and 22 TeV, respectively. 
The integral flux is 1.3±0.4 Crab units at energies above 1 TeV, which is 3.2±1.0 times the 
flux determined by the HESS experiment, i.e., 0.40 Crab unit. 



5. Discussion 

The integrated energy flux above 1 TeV measured by ARGO-YBJ is ~ 1.9 x 10~ 10 
erg cm -2 s _1 , corresponding to a source luminosity, assuming isotropic emission, of L(>1 
TeV)~ 2.3 x 10 34 (D/1 kpc) 2 erg s -1 , where D is the distance to the source. However, due to 
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the limitations in the angular resolution, this flux may also include other contributions apart 
from HESS J1841— 055. Diffuse 7-rays, produced by cosmic rays interacting with matter in 
the Galaxy plane, are expected to contribute to the ARGO-YBJ result. According t o the 



meas urement of diffuse 7-ray flux from the inner Galactic plane using ARGO-YBJ data (IMa 



20111 ). this contribution to the flux from HESS J1841— 055 in the five intervals is less than 
4%. Assuming the HESS shape for the source instead of the symmetrical two-dimensional 
Gaussian shape, the flux would only vary of 2.2%. HESS J1837— 069 is the nearest VHE 
7-ray source with an angular distance of 1.62 degrees. The flux from HESS J 1837— 069 at 



energi es above 1 TeV is 17% that of the Crab with spectral index —2.27 ( lAharonian et al. 



2006b|). Its contributions to the five intervals are estimated to be 5.8%, 2.7%, 1.0%, 0.7%, 
and 0.2%, respectively. The second nearest source is HESS J1843-033, whose flux is still 
unknown. A hot sp ot with a marginal significance of 4.1 s.d. is observed near its position 
( IHoppe et al.ll2007l ). With an angular separation of 2.6 degrees, its contribution is estimated 
to be lower compared with HESS J1837— 069. The contribution from other known VHE 
7-ray sources is negligible . Moreove r , an estimate of the systematic error of ARGO-YBJ 
is described in ( iBartoli et al.l l2012b[ ). With an incomplete list of possible causes, such as 
time resolution variation, event rate variation with environment parameters, and pointing 
error, the systematic error for point sources is found to be less than 30%, and is lower for 
extended sources. Thus, the systematic error of ARGO-YBJ alone is not enough to explain 
the discrepancy. 

HESS J1841-055 is observed by ARGO-YBJ only at high zenith angles (6 > 35.7°), 
while the observation of the Crab is possible also at low zenith angles. This difference may 
cause some systematic errors when estimating the spectrum of HESS J1841-055. To check a 
possible systematic error, observations of the Crab at zenith angles higher than 30° are used. 
With this selection, the average zenith angle is about the same as that of HESS J1841-055. 
The result is that the Crab spectral index varies from (2.58±0.07) to (2.52±0.21), and the 
flux above 1 TeV is (35±28)% higher. Due to the large statistical error, we cannot exclude 
a systematic effect causing the difference of flux. However, even taking this systematic error 
into account, the flux of HESS J1841-055 observed by ARGO-YBJ is still about twice that 
determined by HESS. 

On the other hand, the discrepancy is similar to that found for the two extended sources 
MGRO J1908+06 and MGRO J2031+41 teartoli et alll2012cl lbh. The fluxes measured by 
the EAS arrays Milagro and ARGO-YBJ are much higher than that determined by the 
Cherenkov arrays. Since a good agreement has been achieved on the "standard candle" 
Crab Nebula, some systematic diffe rences between the two techniques should exist only for 
extended sources. As pointed out by lAbdo et al.l (120121 ) . due to their limited FOV, Cherenkov 
telescopes might count the extended emission as background, especially when using the 
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"wobble mode" to estimate this latter. It is worth noting that the "wobble mode" was used 



when HESS observed HESS J1841— 055, and the source was offset by 0.7° ( lAharonian et al. 



20081 ). The "reffected-region technique" is used to estimate the background for spectra in 
a region partially overlapping the extended source. As a result, HESS would measure an 
emission fainter than that measured by ARGO-YBJ. 

Different scenarios have been proposed to explain the emission mechanism of TeV pho- 
tons. VHE 7-rays can be produced via inverse-Compton of background photon fields by 
high energy electrons, or, in hadronic models, by inelastic proton-proton or proton-photon 
interactions. In both scenarios, X-ray and radio synchrotron emissions are expected, there- 
fore the lack of a low energy counterpa rt for HESS J1841— 055 poses the question about 
the nature of the emission mechanism. lAharonian et al.l ( 120081 ) searched for counterparts 
responsible of the VHE 7-ray emission and discussed the possible association with six can- 
didates, marked in Figure 4, which is a zoom of Figure 1 around HESS J1841— 055. Three 
of them are the pulsars PSR J1838-0549, PSR J1841-0524 and PSR J1837-0604, of which 
only the last has a high enough spin-down flux (E/D 2 = 5.2 x 10 34 erg s _1 kpc~ 2 ) to be a 
counterpart candidate. This source is at the boundary of the HESS region, but not far from 



the c enter of gravity of the ARGO-YBJ signal. However, as pointed out in ( lAharonian et al. 



20081 ). since the TeV emission is usually attributed to a relic population of electrons, some 
contribution can be expected also from the other pulsars if they had a much higher spin- 
down luminosity in the past. No catalogued PWNs at longer w avelengths are asso ciated to 
these pulsars, however, according to the recent calculations of ( ITibolla et al.l 120121 ). during 
their evolution ancient PWNs (>>10 kyr) might appear as GeV-TeV 7-ray sources without 
X-ray counterpart. The three other catalogued objects located inside the HESS uncertainty 
region are the SNR G027.4 (also known as Kes73), the high-mass XB AX J1841. 0-0536 
and the diffuse source G26. 6-0.1. The SNR Kes73 lies at the edge of the TeV emission re- 
gion. The point-like nature of AX J1841. 0—0536, the only soft 7-ray source detected within 
the HESS J1841— 055 error ellipse, its variability and the required luminosity (about 10 36 
erg s -1 according to t he ARGO-YBJ data and assuming a distance of 6.9 kpc as inferred 



in (jSguera et al.l 120091 )) exclude its association with the entire emission from the extended 
HESS source. Also the diffuse source G26.6-01, at only 1.3 kpc and well inside the emission 
region, could be responsible at least for part of the TeV flux. In Figure 4 are also reported 
four GeV 7-ray sources from the Fermi-LAT second source catalogue within the extension of 
HESS J1841-055: 2FGL J1839.3-0558c, 2FGL J1836.8— 0623c, a nd the two diffuse sources 
2FGL J1839.0-0539 and 2FGL J1841.2-0459 c dNolan et al.ll2012h. Moreover, the two event 
clusters at energies above 100 GeV found by ( iNeronov fc Semikoz 1 120 10h are shown. Three 
GeV sources are within the two event clusters, suggesting that they may be also VHE emit- 
ters: 2FGL J1841.2-0459c is coincident with the SNR Kes 73, while 2FGL J1839.3-0558c 
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and 2FGL J1839.0-0539 are spatially associat ed to the PSR J18 38-0549 and the diffuse 
X-ray G26.6— 0.1, respectively. As remarked in ( ITibolla et al.l 120121 ). the recent observation 
by Fermi-LAT of GeV sources not firmly associated to X-ray counterparts suggests that 
VHE unidentified sources can be explained as ancient PWNs. 



An hadronic scenario is proposed in (INeronov fc Semikoz II2012I ). These authors consider 
the extended 7-ray emission produced by high-energy cosmic rays escaping from the source 
and diffusing in the interstellar medium (ISM). The 7-ray emission should result from the 
interaction of these cosmic rays with the ISM particles. Such extended emission regions 
should be visible as VHE 7-ray sources with fluxes of order 10 -11 erg cm -2 s -1 above 100 
GeV. From the analysis of the Fermi-LAT data they suggest the young nearby pulsar 
PSR J1841— 0524 as a possible low energy counterpart of HESS J1841— 055. However, as 
already stated, due to the energy balance, this association is not without problems. The 
proton/nuclei contribution to the extended 7-ray flux should generate a comparable flux of 
TeV neutrinos with a spectrum expected to follow the 7-ray spectrum. Thus, the observation 
of high energy neutrinos from the HESS source could provide a crucial test to this model. 
A search for individual neutrino sources over a large fraction of both the northern and 
southern skies has b een carried out by the IceCube detector in the 40-string configuration 
( lAbbasi et al.l 1201 if ). The large background from atmospheric muons reduces the IceCube 
sensitivity to neutrino sources in the southern sky at TeV energies, thus the derived upper 
limits are not stringent enough to constrain the hadronic scenario. Data from the combined 
operation of IceCub e and AMANDA have been used to scan for sources in the Galactic plane 



( lAbbasi et al.ll2012l ) with a neutrino flux sensitivity of about 10 11 — 10 12 erg cm 2 s 1 at 



TeV energies. However, the surveyed range of Galacic longitude (36° < / < 210°) does not 
include the region where HESS J1841— 055 is located. 

In the case of hadronic scenarios one expects the source extension to be much larger 
than seen by Cherenkov telescopes (up to the degree scale). Therefore, the lower angular 
resolution and the large FOV of ARGO-YBJ allows the collection of photons from a larger 
source area. This could partially explain the discrepancy in flux with the HESS Cherenkov 
telescope. 



Recently, iGiacinti et al.l ( 120121 ) found that diffusion of cosmic rays and electrons around 
point sources is strongly anisotropic and shows filamentary structures, which may cause the 
shift of the centroid position between HESS and ARGO-YBJ. 
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6. Conclusions 

Since November 2007 the ARGO-YBJ experiment is monitoring with high duty cycle 
the northern sky at TeV photon energies. Using data up to July 2012, an excess with 
statistical significance of 5.3 s.d. is detected from the direction of the unidentified source 
HESS J1841— 055. The source location and extension are consistent with those determined 
by HESS, however the measured flux above 1 TeV is about 3 times higher. This discrepancy, 
already found in the observation of other extended sources, could origin from the different 
techniques used in the background estimation for extended sources with ARGO-YBJ and 
HESS data. The extended morphology of HESS J1841— 055 and the presence of several 
sources within the 90% confidence error region suggests contributions from more than one of 
them, but so far no clear counterparts in lower-energy wavebands can be identified. However, 
the possibility of a GeV-TeV 7-ray source without any counterpart can not be excluded. 
Both leptonic and hadronic productions of 7-rays have been proposed, but it is not easy 
to distinguish between the two contributions basing only on the 7-ray data. The current 
upper limits to the neutrino flux from the HESS J1841— 055 region are too high to test the 
hadronic model. Further multiwavelength observations from radio to GeV energies and data 
from neutrino telescopes of suitable sensitivity are needed in order to disentangle between 
the different emission possibilities. 
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Fig. 1. — The significance map around HESS J1841— 055 as observed by the ARGO-YBJ 
experiment. The two ellipses for HESS J1841— 055 and HESS J1 837— 069 indicate thei r 
positions and the 68% and 90% contours of their extension regions (lAharonian et all 120 08) . 
The p osition and possible extension of HESS J1843-33 are marked with ellipse ( iHoppe et al. 
20071 ). The stars mark the lo cation of the GeV 7-ray sources around HESS J1841— 055 in the 
second Fermi-LAT catalog (INolan et al.ll2012l ). The solid line indicates the Galactic plane. 
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Fig. 2. — Distribution of 8 2 for the number of excess events around HESS J1841— 055. The 
filled region outline the best fit to simulated data assuming a symmetrical two-dimensional 
Gaussian shape with a = 0.40°. 
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Fig. 3. — Energy density spectrum of HESS J1841— 055 as measured by the ARGO-YBJ 
experiment: the solid line and shaded area indicate t he differential energy sp ectrum and the 
1 s.d. error region. The spectum measured by HESS ( Aharonian et al.ll2008l ) is also reported 
for comparison. Only statistical errors are shown. 
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Fig. 4. — Zoom of Figure 1 around HESS J1841— 05 5. The squares and the dashed circle 
indicate the position of the candidates reported in (lAharonian et al.1 120081 ). The circles 
indicate the two event clusters found in ( Neronov fc Semikoz 201of ) at energies above 100 
GeV. The ellipses and stars are the same as in Figure 1. The solid line indicates the Galactic 
plane. 



